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Abstract 

Most of the conventional models for opinion dynamics mainly account for a fully local influence, 
where myopic agents decide their actions after they interact with other agents that are adjacent to them. 
For example, in the case of social interactions, this includes family, friends, and other strong social ties. 
The model proposed in this contribution, embodies a global influence as well where, by global, we mean 
that each node also observes a sample of the average behavior of the entire population (in the social 
example, people observe other people on the streets, subway, and other social venues). We consider a 
case where nodes have dichotomous states (examples include elections with two major parties, whether 
or not to adopt a new technology or product, and any yes/no opinion such as in voting on a referendum). 
The dynamics of states on a network with arbitrary degree distribution are studied. For a given initial 
condition, we find the probability to reach consensus on each state and the expected time reach to 
consensus. The effect of an exogenous bias on the average orientation of the system is investigated, to 
model mass media. To do so, we add an external field to the model that favors one of the states over 
the other. This field interferes with the regular decision process of each node and creates a constant 
probability to lean towards one of the states. We solve for the average state of the system as a function 
of time for given initial conditions. Then anti-conformists (stubborn nodes who never revise their states) 
are added to the network, in an effort to circumvent the external bias. We find necessary conditions on 
the number of these defiant nodes required to cancel the effect of the external bias. Our analysis is based 
on a mean field approximation of the agent opinions. 
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I. Introduction 

In many complex network problems originating from the social and economic sciences, the aim is to 
study and understand the collective network behavior under particular microdynamics. Opinion dynamics 
models are a subset of these efforts; agent-based models are studied to elucidate the macro-behavior, 
given a prescribed interaction scheme. In opinion dynamics models each node is endowed with a state, 
which adapts and evolves based on observations of the states of other nodes in the network. The state 
is usually a number, which can be considered continuous [1] or discrete [2]. Depending on the problem, 
nodes may represent, e.g., people (3j, firms (4|, or countries (5J. At each time step, each agent observes 
the states of other agents with which it interacts, and then the agents update their states based on these 
observations. 

In this paper we study the effects of different interaction phenomenon on the behavior of opinion 
dynamics models running over networks. Each agent has a binary state ±1. Examples of cases where 
such a model may be applicable are: elections with two major parties (6}, (7J, deciding whether or not 
to adopt a new technology (8), J9j, whether or not to watch a new movie fT0[ , whether or not to follow 



a new fashion trend [11 1, and any other yes/no proposition (go to war or not, referendums, etc.). 

To model the social factors that effect each agent's state, we initially consider two different sources 
of influence: 

• Local influence: Each agent is influenced by its direct neighbors in the network (e.g., family, close 
friends and relatives, and other strong social ties). 

• Global influence: Each agent is influenced indirectly by the average opinion of the population (e.g., 
by observing people on the streets, classes, subway, cinemas, and other social places). For example, 
one takes a sample of the popularity of a fashion trend or a new technological product (such as 
smart phones or laptops) by observing other people in the society. The existence of a strong tie is 
not necessary in this case. Another major example of the global type of influence is the internet. 

We model both of these sources of influence. At each time step, each node observes the fraction of its 
neighboring nodes who have each of the two states. These fractions constitute the local part of the social 
influence, as will be explained in the next section. Each node also observes the fraction of all the nodes 
who adopt each of the two states. The global part of the social influence will be based on these fractions. 

We find an expression for the evolution of the average density of each state over time (i.e. the fraction 
of nodes who adopt s = +1 and who adopt s = — 1) as a function of the initial conditions. We consider a 
network with arbitrary degree distribution and conduct our analysis under the mean field approximation. 
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We also find the probability that all nodes will end up having the same state (i.e., reach consensus), 
and we calculate the expected times to reach consensus on either of the two states. The conventional 



(or "pure") voter model 1 12 1 is a limiting case of our setup, where the global influence is zero, and our 
results can be seen as generalizing those of [13] for the pure voter model. 

Next, we consider the effect of an exogenous field that influences all agents, biasing them towards 
+1. For example, this could model the effects of mass media. We find an expression for the evolution 
of the expected state over time, as a function of the initial conditions. This expression reveals the rate 
of convergence of all nodes to the state s = +1, depending on the magnitude of the external field. 

Finally, we incorporate and study the effect of anti-conformists: stubborn nodes who fix their opinion 
at —1, opposing the external field. In this case, there is never a consensus; instead, the network reaches an 
equilibrium with a certain fraction of agents taking each side. For this model we determine the fraction 
of stubborn nodes required to sway the majority opinion to —1 for an exogenous field of a particular 
strength. 

A. Related Work 

There is a growing literature on opinion dynamics and related models, and we review the most related 
work here. Different opinion dynamics models have different interaction processes and different updating 
schemes, which are contingent on the specificities of the problem. For example, one can choose random 



blocks of adjacent nodes who enforce their opinion on neighbors under certain conditions [14]-[16|. 
In the so-called voter model, at each timestep, a node is randomly chosen, who then copies the state 
of a randomly chosen neighbor (2J, (7J, (T2), (13), (T7J-|[T9j. Conversely, in the invasion process, the 



randomly chosen node imposes its state on a randomly chosen neighbor [20]. 

Learning and trust have been incorporated in the models, where each node has an estimate of how 
reliable its observations might be |j3j , |21|-|25|. In some models, nodes only interact with neighbors 



whose opinions are close to their own [26 ]-[29 [. In [30], agents have inertia in the sense that the longer 



they have had a particular state, the less probable it gets for them to deviate from it. 



In [31 1 — [ 33 1, the effect of stubborn agents (a.k.a., zealots) are studied. These agents commit to a certain 



state and never vary. In 1 3 1 1 the effect of a single zealot is studied, and it is shown that in one and two 
dimensions the system reaches consensus over the state enforced by the zealot, but this is not true in 
higher dimensions. In [32] the voter model with arbitrary number of zealots is examined on the complete 
graph, and lattice graphs in one and two dimensions, and it is shown that the magnetization has Normal 
distribution whose characteristics only depend on the number of zealots, not on the entire population. The 
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upshot is that a few number of zealots can prevent consensus, regardless of the population size. In (33) , 
the expected values of the average state for arbitrary number of zealots on an arbitrary graph is found 
for the pure voter model, and it is shown that consensus is prevented as long as stubborn nodes exist. 
Also, the steady state of the system is shown to be independent of the initial states of the non-stubborn 
nodes. 



In 1 13 1 the pure voter model on heterogeneous graphs is solved under the mean-field assumption. 
This means that nodes with equal degrees are considered to be indistinguishable in dynamics, and the 
connection probability of each pair of nodes is proportional to their degrees. They find the probability to 
reach consensus on each of the states, and also the expected time to reach consensus. In [34], intrinsic 
"flip" rates are heterogeneous, so that some nodes adopt new opinions more frequently than others. 

B. Paper Organization 

The rest of this paper is organized as follows. The basic model, incorporating both local and global 



influence, is formally defined and studied in SectionnTJ Then, Section III studies the effect of an exogenous 



field. Section IV introduces anti-conformists to the model. Finally, we conclude in Section [V] Throughout 



the paper, theoretical predictions are accompanied by the simulation results. 

II. Dynamics in the Absence of Exogenous Bias 

A. Network Model and Basic Dynamics 

We consider a network of N node^j] The network structure is defined through the neighborhoods of 
each node. For node x, let N x denote the neighbors of node x (i.e., the set of nodes adjacent to x), and 
let z x = \N X \ be the degree of node x (the number of its neighbors). The underlying network is assumed 
to be undirected: for two nodes x and y, we have y G N x if and only if x G N y . The network is also 
assumed to be connected: for any two agents x and y, there is a sequence of nodes xq = x, xi, . . . , X£ = y 
such that Xfc G N Xk _ 1 for each k = 1, . . . , I. We will be particularly interested in random network models, 



such as Erdos-Renyi graphs (35) , (36J, preferential attachment graphs [37], and random recursive trees (as 
in (38), (39), with A k = 1). 

Each node has a binary-valued state s x (t) G { — 1,+1}. Nodes have initial opinions s x (0) at time 
t = 0, and time progresses forwards. At each time step, one node is picked randomly, and it updates its 
state. In general, the update will involve the node either keeping the same state, or flipping to the other 



'Throughout, we use the terms "agent" and "node" interchangeably. 
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state. In order to avoid burdensome notation below, we will omit the dependence on time when it is clear 
from the context. Let us denote the state of the system by a vector s = (s±,S2, ■ ■ ■ ,sn)- The probability 
that a node x flips its state at time t is given by 

P{s x (t + St) = s x (t)\s x (t)} = w x (s(t)), 

The first term, with weight 9 G [0, 1], captures the local influence. It is equal to the fraction of neighbors 
who oppose node x. The second part is the global influence, with weight 1 — 9. It equals the fraction 
of the entire population that disagree with node x. Here 9 is a tradeoff parameter between the influence 
of close social ties and the influence of the society. Node x, upon being selected, flips its state with 
probability w x . Hence, if all nodes disagree with x then x flips its state with probability 1. On the other 
hand, if all nodes agree with x then x does not flip its state with probability 1. 

B. A Conservation Law and Consensus Probabilities 

Now we focus on finding the probability that nodes reach consensus on s = +1 and on s = — 1. Note 
that when the system is not unanimous, at each timestep there are nonzero probabilities to move towards 
the states +1 or —1 since, as long as consensus has not been reached, there is at least one node with 
nonzero flipping probability which, by construction, has nonzero probability for being picked to update 
its state. Hence the only absorbing states of the system are those where either all nodes have +1 or — 1. 

We solve the model in continuous-time approximation, examining the system as 5t — > 0, in which 
case for any quantity f(t) we have f( t+St }~fO _>. f = ®L. Let c x denote the expected value of s x . The 
evolution of c x is given by 

c x = -c x + ^ °y+ \^/-J2 c y )■ (2) 

X y&N* \ ally / 

Let nfc denote the fraction of nodes with degree k, and let z = J2 X z x = ]Cfc ^ n k denote the 
average degree of the network. Also define pk{t) to be the fraction of nodes with degree k which have 
s x (t) = +1- Consider the quantity fi(t) = J2k k n kPk> which can be interpreted as the expected value of 
the quantity kpk for a graph with degree distribution {n^}. 

Under the assumption that the degrees of adjacent nodes are uncorrelated, the probability that two 
nodes are connected is proportional to the product of their degrees; i.e., 

P{y G N x \z y = k} = ^, (3) 
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and P{y E N x } = z y /N. Also recall that n k = P{z y = k}. Then by Bayes' rule, we have 



kz 



P{z y = k\y E N x } 



jjf ■ n k _ kn k 



N 



The assumption that the degrees of adjacent nodes are uncorrelated is justified in a number of scenarios. 



For example, this is the case in the classical Erdos-Renyi random graph model [35], |36|. Social networks 
have also been observed to have preferential attachment properties, and in the well-known scale-free 
network model of Barabasi and Albert [37 ] edges are no longer independent since new edges connect 
to high-degree nodes with higher probability. Nonetheless, in Section |II-D we justify that ([3]) still holds 
true for scale-free graphs. 

We will also make a mean-field assumption and approximate that, for a neighbor y E N x with degree 
z y = k, the expected state c y is equal to 2p k — 1, the average state of all nodes with degree k. Under 



these approximations, we have 



°y ~ z x ^{2 Pk - 1) 

k 



kn k 

z 



E 

k 

z 



2kp k n k 



E 



kn k 



Let m(t) = jfYlx s x(t) denote the average state at time t. Then combining (|2j) and 
following differential equation for m(t): 

'2/i 



m{t) = —m(t) + • 



1 +(l-e)m(t). 



(4) 

gives the 
(5) 



Multiplying Q by z x and summing over all nodes gives the following expression for the evolution of 



M(t): 



Compaiing ([5]) with (|6]), observe that 

m(t) 



-(1 



+ 8p(t) = 0. 



Thus, we arrive at the following conservation law: 



^(t) =m(t) ^(1-0) + Op(t) = const. 



(6) 



(7) 



(8) 



For the pure voter model (i.e., with only local influences), we have 9 = 1 and ip = \i is conserved, which 



matches the findings of [19|. 
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Now this conservation law can be used to directly determine the probability that the network reaches 
consensus on a particular state, given the initial opinions at each node. Note that ip equals |(1 + 9) for 
the case where s x = +1 for all x, and ip = ^(1 — 6) when s x = — 1 for all x. For brevity let us define: 



P u d 4 f p{ m (oo) = 1} 



Using the fact that 



P d def p{ m(oo) 



V>(0)= lim ^{t)=P u --{l 

t— >oo Z 



(9) 



■1} 



+ (1-P«).— (1-9), 



the probability of reaching consensus on s = +1 is easily found to be 

pu 

Similarly, we find 



i + 



In terms of the initial states, we can rewrite ( fTO] ) as 

1 1 



2Nz 



and ([IT} also can be rewritten as 

pd _ 



E'-(°)-oAFiEW0). 



(10) 



(11) 



(12) 



(13) 



2 2N ^ w 2A^z 

Figure [I] illustrates the probabilities P u and P d for reaching consensus on ±1 as a function of initial 
fraction of +1 nodes for opinion dynamics run on Barabasi-Albert scale free graphs with 400 nodes. The 
initial opinions at each node are drawn i.i.d., with the initial density of state +1 varying between 0.05 
and 0.95, and the local-global tradeoff parameter 6 drawn uniformly between and 1. As expected, the 
probability P u to reach consensus on +1 is higher when more nodes initially have state +1, and this 
probability is linear in the initial density of nodes with +1. 

Next, we can solve for the dynamics of the average state, m(t). Equations ^ and ^ are a standard 



system of linear differential equations (see |40]-[42] for example). The solution for fj,(t) is 



M (t) 



and for m(t) we have 



m(t) 









+ 




"2 


+ 


z 





Mo - ipo - ~(1 - 0) 



m 



:V>0 + • 



(14) 



(15) 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

p , Initial density of state s=+1 



Fig. 1: The expected probability of reaching consensus on the state +1 for different initial conditions for 
the fraction of nodes with s = +1, on a network of 400 nodes, for different values of 9 between zero 
and one (uniform). The results are averaged over 1000 Monte Carlo simulations. The underlying graph 

with m = 2. 



is a Barabasi- Albert scale-free graph [37 



For both m(t) and p{t) there is a steady-state value and there is a transient exponentially-decaying 
part. Figures [2] and [3] are depictions of the theoretical prediction and simulation results for the steady-state 
values. 




0.4 0.5 0.6 

p , Initial density of state +1 

Fig. 2: Steady state values of m(t) as a function as p(0) for different values of p(0) between 0.1,0.9 
(uniform) and the predicted value of equation Q with 6 = 0.7. The underlying graph is a random 
recursive tree (as in [38 1, with = 1) with 200 nodes. Results are averaged over 100 Monte Carlo 
trials. 



Note that for 6 = (i.e., purely global influence, which is equivalent to a complete graph), equation 
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0.4 0.5 0.6 

Initial density of state +1 



Fig. 3: Steady state values of fi(t) as a function as p(0) for different values of p(0) between 0.1,0.9 
(uniform) and the predicted value of equation Q with 9 = 0.7. The underlying graph is a random 
recursive tree (as in [38 1, with Ak = 1) with 200 nodes. Results are averaged over 100 Monte Carlo 
trials. 



([8]) implies that m is conserved. Also for 9 = 1 the conserved quantity is p. Both of these polar cases 
agree with the results in |13j. 



C. Time to Consensus 

Now let us focus on the expected time required for all nodes to end up with the same state, and 
expected consensus time on s = +1 and s = — 1 individually. Recall that is the fraction of nodes with 
degree k; thus Nrik is the number of nodes with degree k. Of the nodes with degree k, if one changes 
state from —1 to +1, then note that pk changes by l/(Nrik). Let us define 5pk = t^T' an ^ define 



del' 



p {Pk -> Pk + Sp k } 



P k = P {Pk -> Pk~ $Pk} 



(16) 



For the first one we have 
P, 



k - n k (l - p k ) 



1 9 (2p 

2 + 2 l T 



1 + 



-m 



- Pk) 



z 



-711 



(17) 



which implies that a degree-k node with state —1 is pickedand it flips its state at this step. Similarly, 



nkPk 



1 9 (2p 

2 ~ 2 IT 



1 



1 



-m 



nkPk 



1 + 



1 



-P 



-in 



(18) 



Let us denote by T(p) the expected consensus time when the densities for population of nodes with 
various degrees are p = (p\,p2,P3, • • •)■ Note that T(p) is the expected time to reach consensus on either 
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state. We will find the expected time to reach consensus on individual states +1 or —1 afterwards. The 
expected consensus time satisfies the following recurrence relation: 



T(p) =St + J2 [ T (p + ^Pkk)Pk + T (P- 5pkk)P k 



+ 



T{p). (19) 



The first term is the time increment due to an update occurring. The second term accounts for the expected 
time to consensus, given that the selected node has changed its state. The variable k in the sum runs 
over all possible degrees in the network. The last term is the expected time to consensus, given that the 
selected node has made no change in its state. 

Let us normalize time units so that 5t = j^. Define the constants 



x def i_g </, 



def „ 

a = 6 



2ip 



(20) 



Since ift is conserved, we will drop the argument t. 



Let us denote 8 by d k T(p) for brevity. In the continuous approximation (as N — > oo), equation ( [19] ) 
becomes 

J> - Pk )d k T(p) + £ ^ [A + ap k ] dlT(p) = -1. 



(21) 



k k 

To simplify this equation we note that the change in p k , the expected value of p k . We find 



d n - P+ 
dt Pk ~ ^ 



Pi 



n k {\ - Pk) => Pk(t) = A + (pfc(0) - A)e" 



-n fc t 



(22) 



which means that the deviation of p k from the central value A decays exponentially, and thus the expected 
density p k rapidly approaches A for all k. It is observed in the simulations that this phase is very short 
compared to the rest of the process (see Table [I]). Also Figure [4] shows samples of the time behavior of 
p k in simulation. 



Using ( [22] ), the first term in ( [21j ) vanishes, and ( |2T| > simplifies to: 

k 

or equivalently, 

Ef 1 



-6 ib 



1 + 



f)^ + ,-(i-tf)] 2 ^r = -i. 



(23) 



(24) 



Note that when every node has s = +1 then the time to consensus is zero. In this case we have m = 1 
and p = z. Similarly, the time to consensus is also zero when all nodes have s = —1, in which case 
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N 


T 




T x 
T 


100 


61 


3 


0.049 


200 


149 


3 


0.020 


300 


221 


5 


0.023 


400 


317 


6 


0.019 


500 


402 


8 


0.020 


600 


445 


8 


0.018 


700 


496 


9 


0.018 


800 


549 


12 


0.022 


900 


638 


14 


0.022 


1000 


755 


15 


0.020 



TABLE I: Comparing the average time to reach consensus (denoted by T), and the short initial phase 
required for nodes of each degree to reach = A(l ± 0.03) (denoted by T\). The value of 6 is chosen 
randomly in the interval of (0.1, 0.4) for each simulation. T\ is averaged over all degrees. The underlying 



graph is a scale-free Barabasi- Albert graph as in |37|, with m = 2. The quantities are averaged over 100 
Monte Carlo simulations for each case. The initial distribution of states is generated randomly for each 
simulation. 



m = —1 and \i = 0. These observations give us the boundary conditions required to solve the differential 
equation: 

'rty = f(i + 0))=o 
Tty = =g{i-e))=o 



(25) 



Let us define the constant 



del' 



N 



(26) 



02<|!> +1 _02' 

where (z 2 ) = z x i s tne average over all nodes of the degree squared. Then equation ( |2"3"] ) can be 

written as 



A 



1-6 V 
— + I 



1 + 



z , d v - 



(2V) 



The expected consensus time is given by the solution of (23) with the boundary conditions (25 1. The 
solution is readily obtained by integrating twice: 



12 




Fig. 4: Simulation results compared to what 
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k=3 
k=4 

-k=5 
k=6 



predicts. Results are averaged over 1500 Monte Carlo 
simulations (with identical initial conditions), for a Barabasi-Albert graph [37] with N = 500 nodes and 
m = 2. The initial condition was p = 0.7, with nodes selected randomly regardless of degrees. The value 
of 9 is 0.7. 



We can use P u and P d to get a shorter expression for T: 



or equivalently, 



T = A[P U In P u + P d In P d ], 



P u ln P u + (1 - P u )ln(l - P u ) 



(29) 



(30) 



Hence, the time to consensus depends on an entropy-like quantity involving the probability of reaching 
consensus on either state. Note that A < 0. Thus, this expression confirms with the intuition that the 
expected time to reach consensus is largest when the probability of reaching consensus on either of the 
states ±1 is equal, given the initial conditions. 

The theoretical prediction for T is validated in simulation in Figure [5] Also Figure [6] is a depiction of 
consensus time as the population size changes, for a Barabasi-Albert network and random recursive tree 
(as defined in [38], with Ak = 1). 

Now let us examine two polar cases, 9 = 1 and 6 = 0. These polar cases were previously studied 



in [ |T9| and p3| . When 9 = 1, we get the conventional voter model as discussed in [13|. From ( pL2] > 

-N 



and ( 26 1 we have 



A 



pu 



— 7 

z 



-N 



(31) 



13 



which gives 



-N- 



1 



In 1 



+ 



This corresponds to Equation (13) of |13|. 
For the other case, 9 = 0, we have: 



which gives 



-N 



A = -N 
P U = P, 



(1 - p) ln(l - p) + phip 



(32) 



(33) 



(34) 



This matches Equation (25) in (19). Thus, the expressions derived in this expression can be viewed as 



generalizing the previous work [13|, [19| to interpolate between the cases where agents experience either 



purely local or purely global influence. 




p Q , initial density of state s=+1 



Fig. 5: The expected time to reach consensus for different initial conditions for the fraction of nodes 
with s = +1, on a network of 100 nodes, for a Barabasi- Albert graph [37] with m = 2, for different 
values of 6. 



The analysis above found an expression for T, the expected time to reach consensus on any state. 
Next, we refine our analysis by characterizing the time to consensus given that consensus is reached on 
a particular state. Denote by T u the time to reach consensus given that consensus is reached on +1, 



and let T d be defined analogously for consensus on —1. Similar to (19 1, the recurrence relation for T v 
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becomes: 

P u (p)T u (p) =P u (p)5t + 



1 " + P k 



P U {P)T U {P) 



+ Y, [ pU (P + hkk)T u {p + 5 Pk k)P+ + P u {p - 5p k k)T u (p- 5 Pk k)P k 



(35) 



Following a similar line of reasoning that led to (23 1, we arrive at the following differential equation: 



E 



1 



2 + ~ Z 



tp\ ( 1 + 6 ip\ rik 



z N 



6k +-(1-9) 



Using ( [26] ), this further simplifies to: 

A 



1-9 ip 
— + I 



1 



2 z > * 



pUrj-VU 



d 2pu T u 



pu 



pu 



(36) 



(37) 



Now, note that the boundary conditions are different from ( |25| ). In the present case, if all nodes have 
s = — 1 then the expected time to reach consensus on s = +1 goes to infinity. So we have: 

T u (V = 1(1 + 9)) =0 

T u 0= 1(1-9)) ^oo 
Using these boundary conditions, the solution to ((37]) is obtained by integrating twice: 

1 + 9 ip s 



(38) 



1+6* _ ± 
T u = A ^_ g ln 



(39) 



Similarly, for the expected time to reach consensus given that consensus is reached on s = —1, we 



get: 



1-6 I j>_ / 1 
T d = A , 2 n ! ln 



i+e _ ± 

2 z 



-9 V 



(40) 
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Figures [7] and [8] compare the theoretical predictions of (39 1 and (40 1 to simulations. 
We note that one can check that the following holds: 



pUr^U _|_ pdrpd rp 



(41) 



This implies that T is the expected value of consensus, on a probability space of two events, consensus 
on s = +1 and consensus on s = — 1 with respective probabilities P u and P d . Hence, one can also 
use this relation to extract T d , once T u and T are obtained, rather than solving differential equation 
congruent to (37i for T d , which is obtained by replacing all P u s by P d s, and T u by T d . 
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Fig. 7: Consensus time conditional upon consensus over s = +1. The underlying graph is a Barabasi- 



Albert [37| with m = 2 with 200 nodes. Results are averaged over 800 simulations. 
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Fig. 8: Consensus time conditional upon consensus over s = —1. The underlying graph is a Barabasi- 



Albert [37| with m = 2 with 200 nodes. Results are averaged over 800 simulations. 
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D. Dynamics on the Barabasi-Albert Graph 



Let us examine in the passing, the validity of Q for the Barabasi-Albert model |37[ . The Barabasi- 
Albert model is defined through a sequential algorithms for constructing a graph. In each step of the 
algorithm, one new node is introduced and is connected to m existing nodes. (Note that these steps 
are different than time steps in the opinion dynamics models studied above.) Let us say that node j is 
introduced at step j. Then the probability that node j forms a link with an existing node i < j is 

P(i^j)=m tU) (42) 

where Zi(j) is the degree of node % just before node j is added to the network. 

This scheme is called preferential attachment to evoke the notion that existing nodes with higher 



degrees are more likely to receive new links. The denominator in (42 1 is equal to jz(j) = 2mj (i.e., 



twice the total number of links in the network when j is added), so we have 



^) = ~2f =► «(j) = rr^- i . (43) 



The probability of j linking to i becomes 



PM^j) = Z J^M P ( i ^ j ) = z ^i. (44) 

V J > 2mt y J> Nz 



This result is equivalent to saying that the ij-ih element of the adjacency matrix is 1 with probability 
z j^i and is zero otherwise. Using this result, we have 

"2/i 



c v = 2^ i^f c y = L m Nn ^ 2pk ~ 1) = Y 

yt^Nx all y k 

which is consistent with Q. 



2 ^ kn k p k - ^2 kn k 



1 

z 



(45) 



III. Dynamics in the Presence of Exogenous Influence 

Now we extend our model by adding an external bias that tries to align all the states towards +1. This 
external field models the effect of mass media. Examples include advertisements on TV, the internet, 
billboards, theatres, etc. 

Each node assigns a weight to personal observations (which themselves are comprised of local and 
global interactions) and a weight to the external bias. This means that there is a constant probability that 
each node flips towards the external field. We denote this weight by 7 E [0, 1] and the strength of the 
external field is denoted by B G (0, 1]. A larger value of 7 corresponds to less exposure to the media 
bias. It is clear that eventual consensus on +1 is inevitable since B > 0. Now we focus on how the 
dynamics of the system change when the field is introduced. 
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We begin by reexamining the conservation law derived in the previous section. The evolution of the 
expected state of a node at site x is 



c x = 7 



+ - E s v + C 1 



+ (1-7)5. 



Summing up this equation for all nodes gives 



m = (-6>7)m + f^-J /i - 6*7 + 5(1 - 7). 



Multiplying ( 46 ) by z x and summing over all nodes yields 



i l = z l{i-e)m+{i-e) 1 n+ z -[{i-e) 1 + {i- 1 )B). 



(46) 



(47) 



(48) 



From (47 1 and (48 1 we see that ip(t) is no longer a conserved quantity; rather, it evolves as follows 

d 



1> 



dt 



so we have 



^(t) = V(0) + (l-7)7T 5^ • 



(49) 



(50) 



Using this, we can express fi in ( [47] ) in terms of to and time, which yields a first order linear equation 
for m(t). We get: 
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to = —7m. H — —ipo + 7(1 ~~ l)Bt + 5(1 — 7). 



Integrating this equation gives the average state over all nodes 

>(0) 



m(t) = 5(1 - i)t + 



+ 



TO -2-^ + . 

z 



-7* 



For /x(t) we plug ( [52] ) into ( [50] ) to get: 



M«) = y 5(l- 7 )i+[^(0) + |(l 



+ w-V>(o)--(i 



-7* 



(51) 



(52) 



(53) 



Equation ( [22] ) in the presence of the external bias becomes 

1 ~~ l\ 7^(0) 



dt 



Pk = PZ ~P k =n k { 



+ 5 



+ 



l Pk + ^{l-l)Bt 



(54) 



Note that (54) is a standard linear differential equation with integration factor ip(t) = exp(nfc7i) (for 



example, see [43], [44]). We get 

b 



Pk{t) 



7 



t + 



a 



7 nkY 



+ Pk(0) + 



a 



n k T 7 



-n k jt 



(55) 



18 



so the expected densities grow linearly with time and the population is destined to conform to the external 
bias. 

Figure [9] compares the theoretical prediction with simulations. Note that expected states cannot exceed 
+1, so after a finite time we have s x = +1 for all x, which means that the whole population complies 
to the state imposed externally. 




5 10 15 20 25 30 35 

Time 



Fig. 9: m(t) with respect to time, for po = -08 (bottom), 0.14 (middle) and 0.21 (top) for 7 = 0.7 and 
B = 0.04. The underlying graph is a Barabsi-Albert graph [37| with 200 nodes. Results are averaged 
over 500 Monte Carlo trials. 



IV. The Effect of Anticonformists 

Now we seek the minimum number of nodes necessary to constantly adopt s = — 1 in order to cancel 
off the external influence. Let us denote the number of stubborn neighbors of a (non-stubborn) node at 
site x by q x , and denote the average number of stubborn neighbors among regular nodes by q. Let the 
total number of stubborn nodes be Q and let fj, and m be defined as above, only among ordinary nodes. 
Also let z be the average degree, only taking into account the ordinary nodes, not the stubborn ones. 

The equation of motion at (non- stubborn) site x becomes 

1 



7 



-c x + 



Zx + 

Let us define the constants 



-q x + z x 



2(i 

z 



1 



+ 



N + Q 



(Nm - Q) } + (1 - 7)5. 



A 



n 



def 1 

— N 

def 1 

— N 



E 



(56) 



(57) 



y - 
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Then using (56 1 we get 

m = —7 
For /j,(t) we have: 

. 2(1 + 9) . 



e + g 

l + q 



/2# 7 A\ 
m + I — — I n — 70 + 



5(1-7) -jq 



1 + 



+ 



1 + 



To cancel off the terms involving the external field, it suffices to have 

5(1-7) 



(1-0)7 -5(1 -7)" 



This is equivalent to 



7<? 



1 + . 



5(1-7), 



(58) 



(59) 



(60) 



(61) 



so the last terms in (58 1 and (59 1 vanish. This rules out the external excitation from the dynamics. 



Combining (58 1 and (59 1, the equation of motion for the average state becomes 



rh + m7 



9 + q en\ 2 

— — + 1 - — + m 7 
l + q z ) 



+ 



1 + 



,-?)_„ 



1 



1 + 9 



7 2 A(l - 0) - 1 + 



(62) 



This is a second order linear differential equation. To analyze the roots, let us define 

def + q 



1 + 



K 2 d ^ f 1 



en 

z 

1 - 



(63) 



Then the roots of ( 62 1 are 



n 



>'2 



.1 + 9, 

-(Ki + K 2 ) - y/(K! - K 2 f + 4K 3 
2 

-(JCl + K 2 ) + ^(^i - ^) 2 + 4K 3 



(64) 



It is clear that r\ is negative. Now let us examine r 2 . If K\K 2 > K%, then r 2 is negative. This means 
that, if the following relationship holds: 



+ , /) (l-^] >tfA(J 



which is equivalent to 



en 



A < 



+ 



(65) 



(66) 
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then T2 will be negative. First, note that 



N 4^ q x + z x ~ N 



< 



1 \ - 

vE 1 = 1 - 



(67) 



So A is less than or equal to unity. Thus, if we show that the right hand side of ( |66| ) is at least one, 



then ( |66| ) holds. To see this, first note that 

n 



1 z xQx ^ 1 

N ^ q x + z x -N*-" 

X X 



%x — £i 



or equivalently, 



< i. 



Now let us re-write the right hand side of (66 1 for convenience. It is equal to: 



1 



en 



+ 



For the first factor we have: 



Since 9 < 1, this leads to: 



For the second factor, we have: 



So we have found that 



n 



< i 



< 



en 



> 1. 



+ 



1 + -=;>!. 



en ' 



A < 1 < 



+ 



(68) 



(69) 



(70) 



(71) 



(72) 



(73) 



(74) 



which means that (66 1 holds. Thus r<i < 0. 



As we have shown, both roots of ( |62| ) are negative, so the homogeneous answer decays, and the 
particular response prevails. We obtain an expression for the particular response by setting the time 
derivatives in ( |62] ) equal to zero. This gives the steady-state solution: 

9(A(l-6)-l + f) 



m(oo) 



£|(i-f)-0A(i 



(75) 



The denominator is positive as we showed above. For the numerator, note that the following holds: 



A < 1 < 



1 



en 



or equivalently, 



A(l 



<"-t) 



(76) 



(77) 
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So the numerator of (75) is negative, so m(oo) is negative. Remember that the exogenous stimulus was 
trying to bias the states towads s = +1 and the anti-conformists were trying to pull the states towards 

s 



-1. Equation ( |75| ) states that the average steady-state opinion of the whole population is negative, 
which is in favor of the stubborn nodes. 

We conclude that the distribution of stubborn nodes controls the equilibrium bias through A and Q, 



given that (60 1 holds. Figure 10 shows the simulation results and convergence to what d75J) predicts 




Time 



Fig. 10: m(t) with respect to time, for different values of p between zero and 0.9 (uniform), the red line 



is the prediction of (75 1. The underlying graph is random with 500 nodes. The results are averaged over 
100 Monte Carlo trials. 



Now let us examine a special case. Assume that only depends on k, which means that nodes with 
the same degree have equal number of stubborn neighbors. Let the fraction of stubborn neighbors for 
nodes of each degree to be the same; i.e., there is some a < 1 that = ak for all k. This means 

1(1-7) 



a 



(1-0)7-73(1-7)' 



For the parameters in ( 57 ) we have 



A 

n 



a+l 

a 
a+l 



(78) 



(79) 



The final average state becomes: 



m oo 



a+l 



e+az+a 2 z(l-e 
(l+az)(l+a) 



< 0, 



(80) 
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which is negative, meaning that the bias leans towards the stubborn nodes. Note that for the anticonformists 
to be able to succeed, the necessary condition is 

IB 

^MTT- (81) 

which means that the population must care less than some threshold for the inflicted external bias. 



V. Summary and Discussion 

We considered the problem of binary states, or opinions, among nodes that are connected through a 
network of arbitrary degree distribution. We modelled the social influence by two distinct parts. Each 
node finds the fraction of nodes who are adjacent to it and disagree with it, and also the fraction of the 
nodes in the whole network that disagree with it. It then takes a convex combination of those fractions, 
to find a probability, and at each timestep, flips its opinion with this probability. We solved for the 
probability of consensus on either state, and the expected time to reach consensus, as a function of initial 
conditions. We also found the time evolution of the average state over the whole population. 

Then we added an external influence to the model, that tries to bias all nodes towards one of the states. 
We solved for the average state of the system in time. 

To cancel off the effect of the exogenous bias, we introduced stubborn nodes. These nodes constantly 
oppose the external bias. We found the equilibrium average state of the system, and the necessary 
conditions so that the stubborn nodes can succeed in defeating the external field. 

Possible extensions of the model include using a network topology with given degree correlations, and 
using it to improve the mean field assumption. Also, it will be more realistic, albeit analytically somewhat 
formidable, if the influence between each pair of nodes was a function of the distance between them. 
Alternatively, each node, when accounting for the state of other nodes, can assign different weights to 
them, that can be functions of other types of centralities. In this case, the more central a node is, the more 
influence it exerts on the nodes it interacts with. Upon existence of the external field, nodes can be given 
different levels of exposure, or different degrees of resistance against it. This means that, stubbornness 
can be a quality that exists in all nodes, but to various degrees. 

Another useful modification, especially for marketing purposes, is as follows. Suppose s = — 1 means 
that the node have not bought some product, or have not watched some movie. Then, one can modify 
the model so that nodes could not go back once they adopted s = +1 (so-called "susceptible-infected" 
models in the epidemiology literature). This means that, for example, once they saw the movie, they 
cannot un-see it. This will render the model readily testable against empirical data. 
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